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ABSTRACT 



The analysis of sound propagating by multiple paths in 
an ocean at short ranges has been conducted using a Modified 
^ Time Delay Spectrometry (TDS) technique. In this version of 
TDS , a source driven by a linear FM slide and an HP3561A 
Dynamic Signal Analyzer are used to measure the amplitude as 
a function of frequency of signals traveling by different 
paths and having different arrival times. Two sets of data 
from the acoustic test ranges at the Naval Undersea Weapons 
Engineering Station were analyzed for different environmen- 
tal conditions to determine the relative amplitudes of the 
directly propagating and surface reflected signals. Compar- 
isons with simple rough surface scattering theory showed 
reasonable agreement. Results and control software are 
presented and discussed. Recommendations for future 
applications are made. 
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THESIS DISCLAIMER 



The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases 
of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. 



INTRODUCTION 



A . BACKGROUND 

The Naval Undersea Warfare Engineering Station (NUWES) 
in Keyport, Washington, perforins measurements on torpedo 
radiated noise as part of the test and proofing process of 
torpedoes. From these measurements torpedo radiated noise 
baselines are established, torpedoes are accepted or 
rejected on the basis of meeting established radiated noise 
criteria, and design improvements in noise reduction can be 
evaluated. For these reasons, it is important that sound 
transmission loss between the source, the torpedo, and the 
receiver, a system called the Noise Recording System (NRS) , 
be accurately modeled. The current model, based on 
empirical data, appears to be satisfactory under normal 
conditions. However, under unusually rough surface 
conditions scattering effects may disturb the measurement. 

A limited number of experiments and theoretical studies 
have been conducted by Naval Postgraduate School (NPS) 
Students and Faculty to determine the significance of rough 
surface reflection in the multipath acoustic transmission 
loss problem. Brekke [Ref. 1] developed a computer- 
controlled FFT-based dynamic signal analyzer variation of a 
technique called Time Delay Spectrometry (TDS ) to measure 
the separate contributions of multipath propagation. 
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TDS , originally developed by Richard C. Heyser in 1967 
[Ref. 2 ] j- utilizes a transmitted Linear Frequency Modulated 
(LFM) pulse or "swept tone" with constant amplitude. A 
frequency tracking spectrum analyzer uses the difference in 
arrival times to discriminate between the acoustic signal 
traveling by a direct path and those signals traveling by 
reflected paths. 

In a thesis at NPS, LT Ward [Ref. 3] outlined the multi- 
path problem including environmental considerations for the 
NUWES test ranges at Dabob Bay and Nanoose. Further, he 
generated a computer model for predicting transmission loss 
under varying sea states. This model has not yet been veri- 
fied experimentally. However, TDS may now make this 
possible . 

B. OBJECTIVES 

The objective of this thesis is to apply the modified 
TDS multipath measurement system to measure the effect of 
surface roughness on propagation loss as a function of 
frequency and compare the results with theory. The problem, 
theory, and special considerations will be presented first, 
followed by a description of the TDS systems. Analysis of 
data taken at the test range at Nanoose will then be 
presented including a brief discussion of the variability of 
the results. Results, conclusions and suggestions for 
future research will also be discussed. 
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II. PROBLEM DESCRIPTION AND CONSIDERATIONS 



A. THE MULTIPATH PROBLEM 

For the measurement of radiated noise, during part of 
the test run, the torpedo passes by a vertical array of 
three calibrated omnidirectional hydrophones. Maximum 
horizontal range between the torpedo and array is less than 
1000 yards. Torpedo and hydrophone arrays are at 
approximately mid water depths. Water depth is typically 
about 600 feet in Dabob Bay and about 1200 feet at the 
Nanoose range. Transmission loss due to absorption can be 
neglected since geometrical spreading losses are between 40 
and 60 dB greater than absorption at these ranges at sound 
frequencies of interest between 50 Hz and 30 kHz. 

Since the tests are conducted in relatively shallow 
water, acoustic reflection from the boundaries can make 
significant contributions to the total acoustic intensity at 
the receiver. Therefore, in order to model the multipath 
problem adequately, geometry, boundary conditions, and 
environmental effects must be considered. 

1 . Assumptions 

It is assumed that the acoustic signal is of small 
amplitude and that sound propagation can be described by the 
linear inhomogeneous wave equation. Thus the water channel 
acts as a linear filter. It is also assumed that the water 
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column containing the source and the receiver is 
homogeneous. The bottom is considered to be flat and 
horizontal . 

2 . Multipath Geometry 

Figure 2.1 depicts simple multipath geometry. 
According to the approach developed in Albers [Ref. 4:pp. 
49-51] and presented by Ward [Ref. 3:pp. 87-88] and Brekke 
[Ref. l:pp. 25-27], the acoustic signal at the receiver is 
the sum of source signals traveling the direct path and the 
contributions of "image" sources traveling various reflected 
paths. 

Following Brekke' s notation [Ref. l:p. 27], the 
source depth is ZS and the receiver depth is ZR. The 
horizontal separation between source and receiver is R and 
the channel depth is H. In terms of these quantities, the 
direct path distance (XR) , the surface reflected path 
distance (XSR) and bottom reflected path distance (XBR) are 
given by: 



XR = [R 2 + 


(ZR - ZS) 2 ] V 2 , 


(2-1) 


XSR =[R 2 + 


(ZR + ZS) 2 ]V 2 , 


(2-2) 



and 



XBR = [R 2 + ( 2H - ZR - ZS) 2 ] 1 / 2 . 



(2-3) 
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HORIZONTAL RANGE 



SURFACE 




SOURCE XBR - BOTTOM REFLECTED PATH 



Figure 2.1 Simple Multipath Geometry 
(Extracted from Ref. 1] 
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For periodic signals, the contributions of the image 
sources must be considered with respect to both relative 
phase and amplitude. The phases of the reflected signals 
are shifted with respect to the directly propagating signal 
by an amount proportional to the difference in propagation 
path. Additionally the surface reflected signal undergoes a 
complete phase reversal at reflection. The delayed arrival 
times, for a given sound velocity c, are computed as 
follows : 



Surface time delay is 



fc s - 



XSR - XR 



(2-4) 



Bottom time delay is tj-, = 



_ XBR - XR 



(2-5) 



and the time delay between arrivals of surface and bottom 
reflected sounds, or inter-reflected time delay t r , is 



XSR - XBR /0 „ 

t = . (2-6) 

r c 



Assuming an omnidirectional source radiating simple 
harmonic waves, the acoustic pressure wave function at a 
range r is given by: 



P = (A/r) exp j (wt-kr) 



(2-7) 



where A is the pressure amplitude at unit distance, k is the 
wave number (co /c) , and is the angular frequency [Ref. 5: 
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p. 112]. The acoustic pressure at the receiver due to the 
signal arriving via the direct path is given by: 



P d = (A/XR) exp j [ut- (k{ XR} ) ] 

The contribution of the surface "image" is given by: 



(2-8) 



p s = " R s (A/XSR) exp j [out-(k{XSR} ) ] 



(2-9) 



where Rs is the coefficient of reflection for the surface 
[Ref. 3:p. 90]. Similarly, the contribution due to bottom 
reflection is given by: 



where R^ is the coefficient of reflection for the bottom. 

It is possible that Rj., may have an additional phase term 
determined by the nature of the water-bottom composition and 
the angle of incidence [Ref. 3:p. 90]. Assuming that the 
direct and reflected waves are coherent one finds the total 
acoustic pressure at the receiver is the sum of several 
contributions : 



Pb = R b (A/XSR) expj [ cot- (k{ XBR) ) ] 



( 2 - 10 ) 



p tot — p d + p s + Pb * 



( 2 - 11 ) 



Equation (2-11) can be rewritten as: 



p tot = M P d 



( 2 - 12 ) 
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where 



M = 11 - XSR R s ex P(5“t s ) + XBR R b exp(j “ t b )] • < 2 - 13 > 

Ward [Ref. 3:p. 91] defines this factor as the "Multipath 
effect. " 

If the phases of the signals are incoherent, the 
intensity of the combined signals is the sum of the 
intensities of the separate signals, and is given by: 

Itot = Id Is + Ib • (2—14) 

Equation (2-14) can then be rewritten as: 



Itot “ M I Id 



(2-15) 



where 



M Z = [1+I s /I d + I b /I d ] . (2-16) 

Equation (2-16) is defined as the Intensity Multipath 
coefficient. The surface reflection intensity, I s , is 
determined by the condition of the surface. 

The reflection coefficients, R s and R b , are 
determined by the physical properties of the boundaries. R b 
is small at both ranges [Refs. 1, 3] and should be 
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considered significant only when either source or receiver 
is located near the bottom. For purposes of this research 
R b is assumed to be negligible. The surface is assumed to 
provide no mechanism for absorbing acoustic energy and, 
since the specific acoustic impedance of air is small 
compared to that of seawater [Ref 5:p. 127], the energy in 
an acoustic signal incident at the surface is assumed to be 
perfectly reflected or scattered back into the water. 

3 . Surface Reflection 

Urick [Ref. 6:p. 128] describes the sea surface as 
"both a reflector and a scatterer of sound and has a 
profound effect on propagation . . . where source and 

receiver lie at shallow depth." Surface "roughness," 
defined in terms of signal wavelength, determines the nature 
of the acoustic signal after interaction with the surface. 
Urick defines the Rayleigh parameter . kH sin 9, as the 
criterion for acoustic "roughness," where k is the wave 
number (w/c) , H is the rms wave height (peak to trough) and 
is the grazing angle [Ref. 6:p. 129]. When kH sin 9 >> 1, 
the reflection process is called scattering and the sea 
surface randomly scatters all the energy from an incident 
acoustic signal. For kH sin 0 << 1, the surface behaves as 

an acoustic mirror and the reflected wave is completely 
coherent with the incident wave [Refs. 3, 5, 6]. 

The vertical elevation of the sea surface is a time- 
variant random process and is usually modeled statistically. 
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A relatively good first approximation of the probability 
distribution of the surface wave height [Ref. 3:p. 66], 
particularly when kH sin 6 <1 [Ref. 7:p. 344], is the 

stationary Gaussian process. According to developments of 
Beckmann and Spizzichino [Ref. 8:pp. 73-74], the mean time 
average scattering coefficient <q> for a normally 
distributed surface is given by: 



where q 0 is unity (smooth surface), X(k z ) is the spatial 
Fourier transform of the surface probability density 
function or "characteristic function," and k z is the wave 
propagation vector for the acoustic signal. The 
characteristic function for a Gaussian distribution is given 
by: 



where a is the mean-to-peak rms wave height. The 
propagation vector, k z , for the case where the reflected 
angle <P r is equal to the incident angle is given by: 



<q> = X(k z ) q 0 



(2-17) 



X(k z ) = exp (- a 2 k z 2 /2 ) 



(2-18) 



k z = 2k cos 



(2-19) 
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The terra k z is equivalent to the Rayleigh parameter . 

Beckmann and Spizzichino define it as g ^/ 2 and use it as a 
measure of surface roughness [Ref. 8 :p. 82]. Using the 
definition of g and Equation (2-15) , the mean time average 
scattering coefficient of Equation (2-13) becomes: 

<q> = exp(-g/ 2 ) . ( 2-20 

The power reflection coefficient, R^ , is the ratio 
of acoustic energy scattered by a rough surface to the 
energy scattered by a smooth surface (unity) [Ref. 3:p. 68 ] 
and is defined by the expression: 

R = <qq*> = exp(-4k 2 a 2 cos 2 $•!) . (2-21) 

IT 

The surface amplitude reflection coefficient, R s , is then 
given by [Ref. 3:p. 63]: 

R s = R ^/ 2 = exp(-2k 2 0 2 cos 2 $j_) . (2-22) 

Based on experimental determinations, Clay and Medwin [Ref. 

7:p. 344], show that for values of g-*-/ 2 /2 greater than 
unity, Equation (2-19) tends to under-estimate the true 
value of R s . They cite shadowing effects and variation in 
actual wave shapes as reasons for these differences. 
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If it assumed that there is negligible bottom 



reflection, the multipath term of equation (2-12) becomes: 



XR 



M s = [1 ' XSR R s exp(-j W t s )] . 



(2-23) 



The sound pressure level at the receiver due to sound 
travelling by the direct path is given by: 



SPL D = SPLm - 10 log 10 | M s | 2 , 



(2-24) 



where 



SPLm - 20 l°gi0 l p totl • ( 2 - 25 ) 

Equation (2-22) is the total sound pressure level at the 
receiver. 

4 . Environmental Considerations 

Brekke [Ref. 1] and Ward [Ref. 3] include extensive 
environmental data for both NUWES test ranges. The key 
environmental variable affecting multipath propagation is 
extreme variation in the sound velocity profile (SVP) . 

There is the possibility of strong negative gradients 
occurring above 100 ft in depth. The resulting refraction 
may reduce delay time between the signal traveling the 
direct path and reflected paths so that they may not be 
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discernible [Ref. l:pp. 27-37] or may affect the actual 
angle of incidence at the surface. 

B. TDS CONSIDERATIONS 

An analytical discussion of the modified TDS technique 
is presented by Brekke [Ref. l:pp. 16-24]. The points of 
the verification to be considered here are the linear sweep 
or LFM pulse and the response of the HP3561A Dynamic Signal 
Analyzer . 

1 . The LFM Pulse 

Figure 2.2 represents the LFM pulse in the time 
domain. The amplitude of the signal is constant and for 
convenience is assigned a value of unity. The frequency 
range F is the difference between the start frequency f]^ and 
the stop frequency f 2 . The carrier frequency f c is defined 
as the average of f 2 and f^. The sweep has a period of T 
seconds before it repeats. In terms of these parameters, 
the instantaneous voltage of the LFM pulse can be determined 
as a function of time in the expression [Ref. l:p. 17]: 

X(t) = cos[ (tt (F/T) t 2 - TTFt + 2 TTf c t] . (2-26) 

The phase of the signal is: 

4>(t ) = tt (F/T) t 2 + 2 TTfjt . (2-27) 
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Figure 2.2: The LFM Pulse 

(Extracted from Ref. 1) 

The instantaneous frequency f^ is obtained by taking the 
derivative with respect to time of Equation (2-27) and 
dividing by 2tt : 

fj. = (F/T)t + f ± . (2-28) 

The sweep rate (S) is defined as the derivative with respect 
to time of Equation (2-28): 

S= F/T . (2-29) 

In the context of multipath propagation, the 
instantaneous frequency of the signal propagating by the 
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direct path differs in frequency from the reflected signals 
by an amount equal to the sweep rate, S, times the delayed 
arrival time of the reflected signal given in Equations (2- 
4) and (2-5) . 

2 . Response of the HP3561A Dynamic Signal Analyzer 

In the current experimental configuration, the LFM 
pulse is applied to a broad band-width projector. The 
acoustic wave propagates by one or more paths through the 
water to a hydrophone. The electrical signal from the 
hydrophone is recorded in analog form on magnetic tape. The 
recorded signal is subsequently analyzed using the HP3561 
Dynamic Signal Analyzer. The measurements are based on a 
1024 point FFT algorithm. The displayed spectrum contains 
400 bins. Frequency resolution of the analyzer is 
determined by the selected frequency span F of the analyzer 
in terms of the bin separation b: 

b = F/400 Hz . (2-30) 

For a frequency span of 5 kHz, the corresponding frequency 
resolution is 12.5 Hz. The time record length, the 
observation time tQj^ during which the analyzer receives 
data from the signal, is given by: 

t Q bs = Vb sec . (2-31) 
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This is the minimum amount of time required by the analyzer 
to obtain sufficient samples to carry out the FFT algorithm. 
The energy contained in each of the 400 bins is calculated 
and made available for display by the analyzer during the 
observation time. In most operating modes, the HP3561A 
displays only the most recent time record data. In Time 
Capture mode, however, the analyzer acquires and stores data 
contained in 40 contiguous time records. Time Capture will 
be discussed in more detail in the next chapter. 

The consequence of these relations is that the delayed 
arrival times given by Equations (2-4) and (2-5) must be 
greater than t^g in order to resolve the directly 
propagating and reflected signals [Ref. l:p. 48]. For delay 
arrival times less than t^g , the difference between the 
instantaneous frequencies of the direct and reflected 
signals is less than the resolution of the analyzer and 
thus, the energy of both direct and reflected signals are 
contained in the same bins. 

The possibility of energy leaking into an adjacent 
bin and creating additional measurement errors is reduced by 
using a "window function" feature on the HP3561A. The 
result of using a window function is a minimum bandwidth 
separation, f_band, given by [Ref. l:p. 48]: 

f_band = C x F , (2-32) 
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where C, a constant determined by the window function 
selected, is expressed in terms of the frequency span. For 
a Hanning window, C is 0.00375. Thus, if the selected 
frequency is 5 kHz, the resulting f_band is 18.75 Hz. So, 
in order to resolve the direct and reflected signals, it is 
necessary to insure that the delayed arrival times are 
sufficient to maintain the minimum bandwidth separation. 

This is accomplished by setting a lower bound on the sweep 
rate, S m , given by [Ref. l:p. 49]: 

Sm > f_band/t_delay , (2-33) 

where t_delay is given by Equations (2-4) through (2-6) . 

For surface reflection, t_delay = t s . 
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Ill . TPS MULTIPATH MEASUREMENT SYSTEM 



The TDS Multipath Measurement System is comprised of two 
major subsystems [Ref. l:p. 39] : 

(1) Range and test configuration 

(2) Data processing system. 

The range and test configuration described here is specific 
to the data used in this research. The test format was 
developed by John G. Burwell of the Acoustics Division, 
NUWES, Keyport, and has been implemented by Keyport 
personnel in previous TDS research. The data processing 
system presented by Brekke [Ref. l:pp. 42-61] has been 
modified. A description of the hardware configuration is 
provided here. The software used is a slightly modified 
version of the TDS software documented by Brekke [Ref. 1] . 
The key features of the TDS programs will be discussed here 
and software modifications will be documented. 

A. RANGE AND TEST CONFIGURATION 

The data used in this study were obtained at the NUWES 
range at Nanoose, Canada, in February 1986. The acoustic 
projector for the test was the Sonar Acoustic Test System 
(SATS) . The receiver was the Noise Recording System (NRS) 
with three calibrated hydrophones at depths of 330 ft., 390 
ft., and 430 ft. Horizontal separation between SATS and NRS 
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was approximately 300 yds. Wave heights were estimated at 1 
to 2 ft peak-to-peak at the time of the measurements. 

The test signal was a linear sweep from 50 Hz to 30 kHz. 
Three tests were conducted with sweep periods of 1 second, 2 
seconds and 4 seconds. Thirty sweeps were transmitted for 
each test. 

The sound signals arriving at each of the NRS 
hydrophones were recorded on separate channels of a magnetic 
tape recorder. 

B. THE DATA PROCESSING SYSTEM 

1 . Hardware 

The equipment configuration used is shown in Figure 
3.1a. It is essentially the same as that used by Brekke 
[Ref. l:p. 42-43], with an HP thinkjet printer replacing the 
HP2671G printer. For the routine tests, the graphics and 
printing requirements could be satisfied using the HP 
thinkjet as shown in Figure 3.1b. 

The measurements of the signal are carried out by 
the HP3561A Dynamic Signal Analyzer on the recorded signal 
under computer control. Previously, triggering the Dynamic 
Signal Analyzer was done using a 1 pulse-per-second signal 
derived from the output of the time code generator [Ref. 1] . 
This was not possible, since the period of the LFM pulses 
for two of the tests exceed 1 second. To accommodate the 
longer sweep periods, the trigger signal of an HP3314A 
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Arbitrary Function Generator was set to 4 seconds and 
recorded on the edge track of the data tape. 

2 . Software 

Brekke [Ref. 1] provides a detailed documentation of 
the TDS1 program. Here, a general discussion of the 
capabilities of the TDS software will be presented with 
emphasis placed on the features used in this research. The 
modifications used in the current version 3C of the TDS_1 
program will also be documented. 

The original TDS program, TDS_1 , was designed to 
control the HP3561A Dynamic Signal Analyzer in performing 
TDS Multipath Measurements. It is written in HP Basic 
program language 2.0 with extension AP2_1 . The program is 
structured in a modular design illustrated in Figure 3.2. 

It was designed to analyze a series of LFM pulses sweeping 
from 500 HZ to 20 kHz with a repetition period of 1 second. 
The Setanalyzer subprogram in version_l set the HP3561A to 
measure an LFM pulse corresponding to these parameters. In 
order to accommodate the sweep rates and sweep periods for 
the new data described in part A of this chapter, the 
program required some small modifications. The set analyzer 
subprogram which services all the major subprograms was 
altered to allow an operator interactively to select the 
frequency span, start frequency and trigger method of the 
analyzer prior to starting any analysis. 
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The signal is analyzed by one of three methods: 

Quick Analysis, Normal Analysis, and Auto_Analysis . 

Referring to Figure 3.2, one sees that these three methods 
form the three major subprograms of TDS_1. These analysis 
methods are discussed in the following paragraphs, 
a. Quick_Analysis 

In this method the HP3561A is operated in the 
"Time Capture Mode" . The analysis is based on data taken 
from 40 time records stored in the time buffer of the 
analyzer memory. Each time record consists of the data 
processed during one observation time t 0 k s . The total time 
during which the analyzer processes data is 40 x t 0 k s . 

Figure 3.3 illustrates the storage of the LFM pulse in the 
time buffer. The data are displayed as a 3-dimensional 
magnitude map. Each map consists of 60 traces. The overlap 
(or increment) processing selected determines the time 
increment represented by each subsequent trace. For the 20% 
overlap processing used in the TDS program, each trace 
represents a 0.2 x t 0 k s increment of the delay time. The 
analyzer uses 20% new data taken from the subsequent time 
record for each new trace. 

The purpose of the overlap processing is to 
optimize the useful information displayed by the analyzer. 

In selecting an appropriate overlap, the sweep rate, the 
analyzer frequency span, and the frequency span of the LFM 
pulse must be considered. For example, the 4 second LFM 



32 







0123456789 10 40 

HP3561A time buffer records 



Figure 3.3 Storage of an LFM Pulse in the HP3561A 
Dynamic Signal Analyzer Time Buffer for 
a Typical Case 



33 



pulse sweeps through a 5-kHz frequency span in approximately 
678 ms. The time delay between the direct and reflected 
signal is about 100 ms. For the 20% overlap processing used 
in the TDS program each trace represents 0.2 x 80 ms = 16 
ms. Thus, the analyzer uses about 50 traces to display 
direct and reflected signals. If, however, 30% overlap 
processing is used, each trace then represents 24 ms. In 
this case only 35 traces are used to display direct and 
reflected signals, leaving 25 traces displaying no useful 
data. If 10% overlap processing is selected, each trace 
represents 8 ms. The number of traces required to display 
direct and reflected signals is approximately 90 which 
exceeds the 60 traces that can be displayed by the analyzer. 

Quick_Analysis provides the operator a quick 
overview of a "Spectral Segment" of a single sweep of the 
signal to be analyzed. A "spectral segment" is the portion 
of the LFM pulse that can be displayed by the analyzer. If 
the start frequency of the analyzer is 500 Hz and the 
frequency span selected is 5 kHz, then only the "segment" of 
the LFM pulse between these two frequencies can be 
processed . 

From the spectral 3-D map, the arrival times of 
the "Spectral Segment" for both direct and reflected signals 
can be determined by adjusting the start_time of the 
analyzer. The time when the reflected signal vanishes from 
the "Spectral Segment" should also be determined. These 
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values are used as input parameters for NormalAnalysis and 
Auto_Analysis . After the determination of the arrival times 
and "upper limits of the delayed start_time" , or time when 
the reflected signal vanishes, the program displays 
individually all magnitude maps between these time limits 
for operator input. The operator selects peaks for direct 
and reflected signal using the cursor on the Spectrum 
Analyzer. These values are recorded by the program and are 
printed out at completion of the analysis. 

It should be noted here that Quick_Analysis 
requires only a single trigger since only one sweep is 
captured. The requirement for a sequence of synchronous 
triggers exists when Normal_Analysis or Auto_Analysis will 
be performed. 

b. Normal_Analysis 

In this method 8 sweeps are RMS averaged to 
reduce effects of signal fluctuation and increase signal to 
noise ratio before making measurements. The measurements 
are made in a manner similar to those in Quick_Analysis in 
which the operator must select peak values using the cursor 
on the HP3561A. The program increments the delay_step after 
each measurement to measure the frequency bin of the later 
arrival time. These measurements are repeated for every 
delay_step in the time interval determined by the 
Quick_analysis method. 
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The delay_step should increment each measurement 
by a frequency larger than the frequency resolution of the 
analyzer. 

c. Auto analysis 

This method automates the measurements performed 
by Normal_Analysis method. Due to a high noise level at the 
lower end of the spectrum, the operator's observation was 
necessary to discriminate between direct and reflected 
signals. Consequently, the Auto_Analysis method was not 
used in this study. 

d. Analysis Products 

Each of the methods produces a table of voltage 
levels for direct and reflected signals as well as a noise 
level. The program enables the operator to plot the 
magnitude of direct and reflected signals as a function of 
frequency. This requires that all the tabulated frequency 
values be monotonically increasing. This is not always the 
case, since often the reflected path peak is recorded at a 
lower frequency than the previous measurement due to the 
uneven nature of boundary reflection. Using a longer delay 
time will reduce this problem, but with the result that 
fewer samples are recorded. 
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IV. DATA ANALYSIS 



The analysis of a "spectral segment" of LFM pulse with 
the 4 second period is presented here. Procedures used here 
are the same for the analysis of any "spectral segment" of 
any LFM pulse. Multipath problem considerations are 
discussed with respect to TDS parameters. Measurements made 
using the Normal_Analysis method are used to determine the 
effects that surface reflection or scattering have on the 
acoustic field at the receiver. Measurements by the 
Quick_Analysis method on several different sweeps are 
compared to the Normal_Analysis results to estimate the 
sweep to sweep variability in results. 

A. ANALYSIS PARAMETERS 
1 . Time Delays 

The sound velocity profile given in Figure 4.1 is 
typical for the Nanoose range in February [Ref. 9]. Because 
the effects of refraction due to this profile are small for 
horizontal ranges of a few hundred yards, the water column 
will be assumed to have a constant sound speed of c = 4850 
ft/sec. The multipath parameters for the problem are given 
in Table 1 . 

The minimum time delay given by Table 1 is 82 ms. 
The minimum frequency span setting for the HP3561A is 5 kHz 
corresponding to t^g = 80ms. 
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TABLE 1 



DISTANCES AND TIME DELAYS 

t s (ms) 
8 2 ms 
91 ms 
96.4 ms 

Source depth: 300 ft 

Horizontal Range: 100 yds (300 ft) 

c = 4850 ft/sec. 



RCVR (depth) 
upper (330) 
middle (390) 
lower (430) 



XR (yds) 
100 yds 
104 yds 
210 yds 



XSR (yds) 
233 yds 
251 yds 
266 yds 



2 . The Rayleigh Reflection Coefficient 

Assuming transmission loss to be due to only 
# spherical divergence, the sound pressure level at the 
receiver due to the direct signal is given by: 



SPL S = SPL[ P ( 1) ] - 20 log XR . (4-1) 

The sound pressure level due to a surface reflected signal 
is assumed to be given by: 

SPL r = SPL [P(l)] - 20 log XSR + 10 log (R s ) . (4-2) 

P(l) is the sound pressure at unit distance from the source, 
subtracting (4-2) from (4-1) yields: 
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10 log R s = SPL d -SPL r + 20 log (XSR/XR) . 



(4-3) 



If all terms on the right hand side of (4-3) are 
known, then a value for R s can be determined. Using 
equation (2-21) , one can solve for the mean-to-peak wave 
height. For six sweeps analyzed using the Quick_analysis 
method, these calculations gave an average value for o of 
0.54 ft at frequencies less than 1200 Hz. The standard 
deviation at these frequencies was on the order of 0.05 ft. 
At frequencies above 1200 Hz, the calculated values for o 
decreased . 

These results are consistent with data given by Clay 
and Medwin [Ref. 7:p. 344]. Since the wave heights were 
estimated between 1 and 2 ft, 0.54 feet for o seems tcf be 
reasonable. So for the data used here, the measurements made 
by the Modified TDS Measurement System appear to be valid. 

3 . Selecting Analyzer Settings 

Since the minimum time delay is 82 ms (Table 1) , the 
maximum observation time t 0 j-, s which permits the HP3561A to 
resolve the direct and reflected signals is 80 ms. This 
corresponds to a 5-kHz spans on the analyzer. As pointed 
out in Chapter II, choosing the maximum observation time 
gives the best frequency resolution on the analyzer. 

Figure 4.2 is the 3-D map output of Quick_Analysis . 
The start frequency selected here is 500 Hz. Ideally, each 
5-kHz span in the sweep would be examined successively. 
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B. MEASUREMENTS 



The received Acoustic Level RL(f) at the face of the 
hydrophone is given by [Ref, l:p. 74]: 

RL ( f ) = ML(f) -CAL(f) -HS (f) +IL(f ) -47-AG 



where : 



ML(f) is measured level in (dB re IV) 

CAL(f) is measured level of the calibration signal 
HS ( f ) is hydrophone sensitivity 
IL(f) is calibration insertion loss 
AG is the Auto Gain of the Amplifier. 

The levels ML(f) for these measurements are provided 
in outputs from the TDS_1 program in Appendix B. For 
selected frequencies used in this thesis, HS(f) and IL(f) 
are found in Table 2. These values are from the calibration 
records for the equipment provided by NUWES. The values for 
AG and CAL(f) used here are +40dB and OdB. Since kH sino is 
approximately unity at 1.2 kHz, some assumptions about 
signal reflection at the surface can be made. For f < 1.2 
kHz, the effect of specular reflection from the surface is 
determined by Equation (2-22) . The upper and lower bounds 
of the acoustic pressure at the receiver are given by: 
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TABLE 2 



HYDROPHONE SENSITIVITIES AND INSERTION LOSSES* 

Hydrophone 



Selected 


Upper 




Center 




Lower 




Frequency 


HS 


IL 


HS 


IL 


HS 


IL 


1 kHz 


-181.6 


-20.0 


-182 . 0 


-20.1 


-182 . 2 


-20.1 


2 kHz 


-181.0 


-20.0 


-181.5 


-20.1 


-181.0 


-20.1 


3 kHz 


-181.2 


-20.0 


-182 . 1 


-20.1 


-181.7 


-20.2 


4 kHz 


-182.4 


-20.0 


-182.4 


-20.2 


-182 . 3 


-20.2 


5 kHz 


-182.3 


-20.1 


-182 . 0 


-20.2 


-182 . 3 


-20.3 


*Data 


provided 


by NUWES 











E-max = p Dd + xsR R s ) ( 4 “ 4 ) 

and 

VD 

p min = P D d ~ xH R s } • ( 4-5 ) 

For f > 1.2 kHz, the effect of scattering and 

incoherent addition from the surface on the average acoustic 
intensity at the receiver is given by Equation (2-15) where: 



Mi = [i + i s /i d ] 



(4-6) 



levels 



The ratio I s /I r is equal to the ratio of measured 
of both direct and reflected signal. Since the 
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Normal_Analysis method uses data averaged over several 
sweeps, these values should be used to calculate Mj . 

Multipath coefficients determined for the selected 
frequencies in Table 2 are displayed in Table 3 . 



TABLE 3 

MULTIPATH COEFFICIENTS AT SELECTED FREQUENCIES 



f kHz 
1 
2 

3 

4 

5 



type 

coherent (M s ) 
incoherent (Mj) 
incoherent (Mj) 
incoherent (Mj) 
incoherent (Mj) 



Coefficient 
0.58 to 1.42 
1.53 
1.17 

1.17 

1.18 



C. VARIABILITY OF RESULTS 

Figure 4.3 shows a comparison of the results of 
Quick Analysis and Normal_Analysis for the same conditions. 
Normal Analysis, as discussed earlier, makes an RMS average 
of 8 sweeps. In order to check the variability of results 
Quick_Analysis was performed on six different sweeps. The 
energy averaged level of the measured band levels is given 
by: 



<X> = 10 log 10 t £io lvl /10/ N ] , 



(4-7) 
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where LVL is the measured level in a given frequency band 
and N is the number of measurements. Table 4 gives the mean 
level <X>, and "normalized" level of the standard deviation, 
10 log(S) - <X>, for several frequencies. 

The Rayleigh reflection parameter varies from about 
1 to 5.5 over the frequency span in the experiment. The 
energy averaged Quick_Analysis measurements are close to the 
measurements obtained by the Normal_Analysis method. This 
indicates that the Normal Analysis measurements provide a 
fairly good approximation to the average acoustic energy in 
each band of the measured signal. 



TABLE 4 

STATISTICAL COMPARISON OF SURFACE REFLECTED LEVELS 

(dB re IV) 



Freq (Hz) 


Quick Analysis 
<X> 


Normal Analysis 


101og 10 (s) 


1000 


-12 . 6 


-12 . 0 


-13.4 


1500 


-13.9 


-14 . 6 


- 3.8 


2000 


-18.4 


-18.3 


- 0.25 


2500 


-25.7 


-26.1 


- 2.4 


3000 


-21.2 


-21.9 


- 1.5 


4000 


-24.0 


-21.5 


- 4.0 


5000 


-24 . 4 


-22.9 


- 2.7 
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V. DISCUSSION AND CONCLUSIONS 



A. DISCUSSION 

1 . Procedure 

Two factors which caused the greatest difficulty in 
analyzing data were the way in which the trigger was set and 
the choice of the number of sweeps for each data run. As 
discussed earlier, the Dynamic Signal Analyzer was triggered 
by signals recorded on the edge track of the data tape. The 
Honeywell 5600E magnetic tape recorder has a total of 14 
channels. This includes 13 data channels and one edge 
channel, typically used for voice comments. The edge 
channels can be recorded without disturbing the data 
recorded on the data channels. The edge tracks, located on 
the "edge" of the tape, sometimes suffer degradation caused 
by tape guides on the recorder. The trigger for the HP3314 
Arbitrary Function Generator was set to 4 seconds and 
recorded on the edge track. Several attempts were made to 
synchronize the trigger with the start of the LFM pulse. 

The data analyzed in this thesis required approximately 1.9 
second trigger delay. For the overlap processing used in 
the TDS program, only 80 percent of the LFM pulse was 
processed in Quick_Analysis . Since the trigger was not 
synchronized to the start of the LFM pulse, 20 percent of 
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the signal between 8 kHz and 12 kHz could not be captured by 
the analyzer. 

The trigger level remained high for approximately 
1.7 seconds. This made Quick_Analysis possible for parts of 
the sweep arriving during times of high trigger. 
Normal_Analysis was not possible, however, because the 
analyzer continued to trigger and sample different portions 
of the LFM pulse. 

For further TDS measurements a sharp, short-duration 
trigger is necessary and should be recorded on tape with the 
data. The data Brekke used had a one second sweep period 
which was synchronous with the range time signals. A 1-Hz 
trigger was derived from the IRIG-B, time code generator 
translator. If a one second sweep period is maintained, the 
frequency span of the signal will have to be reduced to 
accommodate a slower sweep rate. If a slower sweep is used, 
such as the 4 second sweep used here, the advantage of a 
larger frequency span is gained. In this case, a trigger 
signal synchronous with the start of the LFM pulse should be 
recorded . 

The total data time for all three tests was 3.5 
minutes, with two minutes the maximum total time for slowest 
sweep rates. Normal_Analysis requires 8 sweeps to make one 
measurement. For thirty sweeps, a maximum of 3 measurements 
can be made before the tape must be rewound, repositioned 
and started again. For the 5-kHz span used here 
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approximately 42 measurements are required and consequently 
the rewind process was repeated a minimum of 14 times. With 
five minutes of data recorded, the rewind process would be 
repeated only 5 times. Since the rewind process totals 
about 3 minutes the analysis time would be reduced 
significantly. An additional benefit is that increased 
information available from the longer record could be used 
to determine the variability of the surface reflection. In 
future tests, the data recorded should be of a longer time 
at one sweep rate with 5 minutes of data accumulation as a 
minimum. 

2 . Software 

The modifications made to the TDS_1 program have 
made it more flexible in making analysis of data. 

Eliminating a requirement for the plotter shown in Figure 
3.1 (b) , would streamline the process considerably. 

Although HP Thinkjet graphic printouts are not as detailed 
or colorful as the plots of the HP7470A plotter, for 
perishable measurements the HP Thinkjet products are 
adequate . 

The Quick_Analysis printouts given in Appendix B.l 
read in 4 ms increments; this corresponds to 16 ms of time 
in the Normal_Analysis printouts, Appendix B.2. 

The program in current form can be used for routine 
measurements. Minor modifications need to be implemented 
for the changes just discussed, but personnel who would use 
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the products and conduct the measurements should determine 
the final form of the program. 

B. CONCLUSIONS 

The Modified TDS Measurement System presented by Brekke 
[Ref. 1] provides useful information about an acoustic 
signal propagating by direct and surface reflected paths. 
When the sea surface has some degree of roughness, 
measurements by this method permit the effect of the surface 
reflected signal and the direct signal to be considered 
separately. Further modifications to the measurement system 
are still required. In its current configuration, the 
system is cumbersome. The removal of the requirement for a 
plotter could streamline the analysis for routine tests. 

The HP3561A Dynamic Signal Analyzer can perform an 
energy average of a recorded bandwidth. This technique may 
be more accurate than the current method of recording only 
the peak value, particularly for Quick_Analysis. In 
Normal_Analysis the averaging over 8 sweeps, eliminates most 
of the variability so there should not be a significant 
difference between the average of peaks and the energy 
average. In Quick_Analysis , however, the averaging over the 
bandwidth would improve the accuracy of measurements, 
particularly for the reflected signal at higher frequencies. 
The results obtained here are not conclusive. More data 
needs to be analyzed. The results reported here were made 
for a limited frequency range. Results for different sweep 
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rates and frequency spans should be investigated and 
compared . 

In any future tests, more extensive measurements and 
recording of environmental conditions need to be made. If 
reflected signals are distinguishable and if more than one 
is present, the operator may use TDS_1 to compare any of the 
reflected signals with the directly propagating signal. 
Future modification to the TDS program should allow for the 
analysis of more than one reflected signal. 
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APPENDIX A 



RANGE AND TEST DATA 



This list provides a guideline of information to be 
recorded at the range just prior to conducting tests using 
the Modified TDS Measurement System. 

1. ENVIROMENTAL 

(a) Wind-Speed and Direction 

(b) Sound Velocity Profile 

(c) Wave Heights 

2. CONFIGURATION OF TEST 

(a) Horizontal Ranges Between Source and Receivers 

(b) Depths of Source and Receiver 

(c) Range to Nearest Land 

(d) Water Depth at Test Sight 

3. TEST PARAMETERS 

(a) Calibration Signal 

(b) Auto Gain Settings 

(c) LFM Pulse 

- Frequency Span 

- Repetition Period 

- Number of Sweeps 

4. CALIBRATION RECORDS FOR RECEIVING EQUIPMENT 

(a) Hydrophone Sensitivities 

(b) Insertion Loss Data 
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APPENDIX B 



TABULATED MEASUREMENT RESULTS 
NPS : UPPER HYDROPHONE 



ANALYIS METHOD: QUICK ANALYSIS 

DATE OP ANALYSIS: 12/1/86 



RECORDED MAGNITUDE LEVELS ( dB re IV) 



DELAY( MS ) 


DIRECT MAGN 


DIRECT PRQ(KZ) 


REFL MAGN 


REFL PRQ( HZ ) 


NOISE MAGN 


1900 


-9.33 


500 


19.75 


0 


-51.16 


1904 


-9.21 


563 


19.75 


0 


-49.99 


1900 


-8.85 


662 


19.75 


0 


-49.69 


1912 


-7.76 


750 


19.75 


0 


-49.51 


1916 


- 9.28 


075 


19.75 


0 


-40.67 


1920 


8.09 


1013 


19.75 


0 


-47 . 57 


1924 


-8.46 


1125 


19.75 


0 


-45.92 


1928 


- 6 . 39 


1238 


19.75 


0 


-43.42 


1932 


-5.06 


1300 


-14.91 


550 


-42.00 


1936 


-6.35 


1363 


-14.81 


725 


-42 . 92 


1940 


-9.85 


1512 


-13 . 85 


707 


-44.01 


1944 


-11.53 


1575 


-14.13 


808 


-43.88 


1940 


-9. 96 


1725 


-12.80 


1038 


-42.29 


1952 


-9.93 - 


1800 


-10.42 


1.125 


-40.59 


1956 


-11.09 


1913 


-8.09 


1263 


-39.72 


1960 


-12 . 27 


2050 


-7.14 


1300 


-40.14 


1964 


-11.24 


2113 


-9.77 


1413 


-41.08 


I960 


-12 . 70 


2263 


-10 . 39 


1438 


-42 .77 


1972 


-11.76 


2325 


-17.93 


1563 


-43.72 


1976 


-13.69 


2438 


-14.58 


1780 


-44.94 


1980 


-15 .24 


2513 


-11.66 


1825 


-44.79 


1984 


-15. 29 


2675 


-13 . 25 


1862 


-44 . 20 


1988 


-13.84 


2775 


-15.29 


2050 


-43.45 


1992 


-12.60 


2875 


-15.05 


2150 


-43 . 15 


1996 


-11.92 


2963 


-13.69 


2188 


-43.06 


2000 


-11 . 79 


3075 


-15.71 


2325 


-43.49 


2004 


-11.76 


3138 


-21.23 


2380 


-43.11 


2008 


-14.02 


3250 


-20.03 


2525 


-43.25 


2012 


-14. 14 


3375 


-21.91 


2662 


-43.60 


2016 


-13.78 


3463 


-22.64 


2687 


-43.72 


2020 


- 14.54 


3550 


-19.82 


2863 


-44. 14 


2024 


-16.51 


3650 


-20.56 


2925 


-44.03 


2028 * 


- 16.61 


3775 


-20.28 


3088 


-44.07 


2032 


-16.57 


3063 


-19.28 


3080 


-43 . 96 


2036 


-16.46 


3950 


-21.54 


3175 


-43.60 


2040 


- 17.92 


4050 


-22 . 95 


3325 


-43.76 


2044 


-17.19 


4188 


-20.74 


3513 


-43.06 


2040 


-16.93 


, 4288 


-19. 11 


3525 


-42.96 


2052 


-17.07 


4375 


-25.74 


3625 


-42.80 


2056 


-17.51 


4475 


-18.46 


3788 


-43.04 


2060 


-17.06 


4588 


-19.91 


3875 


-43.03 


2064 


- 17.76 


4675 


-21.03 


3900 


-43.03 


206 8 


-18.31 


476 3 


-20.02 


4113 


-42.63 


2072 


- 18 . 96 


4913 


-21.88 


4150 


-42 . 52 


2076 


-18.96 


5012 


-24.44 


4250 


-42.66 


2080 


-18.49 


5150 


-25.28 


4425 


-42 . 09 


2084 


-16.71 


5262 


-24.66 


4463 


-43.21 


2088 


-15.13 


5313 


-21.20 


4625 


-43.59 


2092 


- 16.02 


5388 


-23.54 


4662 


-44.61 


2096 


-16.85 


5500 


-25.07 


4750 


-44.73 


2100 


-21.43 


5500 


-27.18 


493 8 


-44.39 


2104 


19.75 


0 


-22.31 


5038 


-45.11 


2108 


19.75 


0 


-21.15 


5125 


-45.06 


•% i i *» 


l Q -7S 


O 


-20.63 


5188 


-44.84 
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NPS 



MIDDLE HYDROPHONE 



ANALYIS METHOD: QUICK ANALYSTS 

DATE OF ANALYSIS: 11/30/86 



RECORDED MAGNITUDE LEVELS ( dB re IV) 



DELAY( MS ) 


DIRECT MACN 


DIRECT PRQ( HZ) 


REFL MACN 


REFL PRQ(H7) 


NOISE MAGN 


1930 


— 6 . 60 


563 


19.75 


0 


-47.81 


1934 


-7.00 


575 


19.75 


0 


-46.91 


1938 


-4.90 


750 


19.75 


0 


-46.20 


1942 


-4.90 


775 


19.75 


0 


-45.30 


1946 


-6.52 


925 


19.75 


0 


-43.57 


1950 


-5.85 


1075 


19.75 


0 


-41.11 


1954 


-4.62 


1213 


-10.46 


537 


-38 . 72 


1958 


-2.35 


1275 


-11.74 


613 


-36.62 


1962 


-2. 17 


1350 


-12.51 


775 


-36.36 


1966 


-3.96 


1400 


-8.80 


838 


-37 . 98 


1970 


-7.69 


1575 


-10.68 


912 


-39.35 


1974 


-7.91 


1713 


-9.84 


1075 


-38.47 


1978 


-7.56 


1775 


-7.72 


1138 


-37.48 


1982 


-7.83 


1250 


-7 . 83 


1250 


-37 . 80 


1986 


- 11 . 70 


1950 


-8.40 


1313 


-39.74 


1990 


10.60 


2063 


-12.51 


1363 


-41.76 


1994 


- L2 .05 


2125 


-19. 36 


1525 


-42 . 98 


1998 


- 14.3 1 


2225 


-15.91 


1713 


-42 . 82 


2002 


- 18.56 


2375 


-14. 30 


1800 


-42.52 


2006 


- 17.16 


2563 


-11.77 


1913 


-41.81 


2010 


14.03 


2625 


- 10. 19 


1975 


-40. 39 


2014 


- 11.78 


2712 


- 12 . 24 


2075 


-40. 37 


2018 


-9.89 


2800 


-18.73 


2162 


-41.09 


2022 


-9 . 97 


2888 


-18.77 


2263 


-42 . 49 


2026 


11.22 


2975 


-26 . 56 


2388 


-42 . 87 


2030 


- 12 . 62 


3075 


-21.96 


2513 


-42 . 76 


2034 


14.96 


3250 


-20.09 


2600 


-42 .19 


2038 


-12.72 


3325 


-20.20 


2700 


-41 . 66 


2042 


-12.21 


3438 


-17.33 


2787 


-40. 95 


2046 


-12.87 


3488 


-21 . 88 


2838 


-40.82 


2050 


-14.32 


3625 


-20.61 


3075 


-40.81 


2054 


-13.76 


3738 


-16.74 


3100 


-40 . 84 


2058 


-13.21 


3825 


-18.80 


3238 


-40.11 


206 2 - 


- 14.37 


3938 


-13.67 


3337 


-39.70 


2066 


- 13.54 


4025 


- 16.05 


3400 


-39.76 


2070 


- 12.61 


4150 


-22 . 63 


3475 


-39.78 


2074 


12.24 


4213 


-17.31 


3663 


-39.95 


2078 


-11.85 


4350 


-16.63 


3738 


-38 . 92 


2082 


- 12.54 


4400 


-15.81 


3825 


-38 . 12 


2086 


-12.14 


45 50 


-16.08 


3938 


-37.79 


2090 


-12.81 


4625 


- 17.19 


4088 


-37. 49 


2094 


-12.88 


4750 


- 14.96 


4188 


-37.41 


2098 


- 13.14 


4850 


-13.11 


4250 


-38 . 29 


2102 


-13.76 


4950 


-13.13 


4350 


-38.96 


2106 


- 13.53 


5038 


-13.67 


4400 


-38 . 23 


2110 


-12.81 


5200 


-15.56 


4550 


-38.56 


2114 


-12.17 


5225 


-16.06 


4588 


-39.12 


2118 


-15.67 


5288 


-18.46 


4763 


-40.04 


2122 


-15.49 


5463 


-19.72 


4825 


-40. 48 


2126 


- 17.01 


5500 


-22.65 


5012 


-41.58 


2130 


19.75 


0 


- 20.71 


5012 


-42.47 


2134 


-20.61 


563 


-22.11 


5200 


-42 . 59 



54 



NPS 



LOWER HYDROPHONE 



ANALYTS METHOO: QUTCK ANALYSTS 

OATS OF ANALYSTS: 11/30/86 



RECORDEO MAGNTTODE LEVELS ( dB re IV) 

OELAY(MS) OTRECT KAGN OTRECT FRQ( HZ ) REEL MAGN REFL PRQ( HZ ) NOISE MAGN 



I960 


19.75 




0 


19.75 


0 


-53.79 


1954 


“12.32 




725 


19.75 


0 


-53 . 47 


1958 


*11.40 




800 


19.75 


0 


-52.60 


1962 


-12.43 




925 


19.75 


0 


-51.99 


1966 


“12 . 62 




1063 


19.75 


0 


-50.91 


1970 


-10.87 




1150 


19.75 


0 


-50.04 


1974 


-8.40 




1275 


-19.97 


537 


-48 . 96 


1978 


-7.57 




1325 


-20.77 


675 


-48 . 13 


1982 


-9.42 




1388 


-19.33 


800 


-48 . 85 


1986 


-13 . 17 




1500 


-17.37 


850 


-49.74 


1990 


-14.19 




1688 


-17. 16 


1000 


-49.13 


1994 


-12.49 




1738 


-16.64 


1075 


-48.17 


1998 


-13.37 




1813 


-16.99 


1150 


-48 . 32 


2002 


-15.20 




1925 


-19.58 


1275 


-49.24 


2006 


-15.11 




2050 


-20.34 


1363 


-49. 33 


2010 


-16.04 




2138 


-23.41 


1475 


-49.32 


2014 


-17.56 




2238 


-21 . 40 


1600 


-49.20 


2018 


-18.68 




2375 


-19.95 


1725 


-49.21 


2022 


-19.70 




2425 


-16.20 


1775 


-49.95 


2026 


-19.58 




2600 


-17.47 


1837 


-49.89 


2030 


- 18.23 




2700 


-21.19 


1950 


-49.96 


2034 


-16.74 




2775 


-26. 30 


2025 


-49.72 


2038 


-15.37 




2875 


-25 . 22 


2225 


-49. 12 


2042 


-15.73 




2963 


-27.51 


2325 


-49.24 


2046 


-15.33 




3075 


-30.01 


2425 


-49.31 


2050 


-16 . 37 




3138 


-26.34 


2487 


-48 . 66 


2054 


-17.29 




3300 


-25.28 


2588 


-48 . 14 


2058 


-16.97 




3400 


-24.01 


2725 


-47.54 


2062 


-17.63 




3463 


-21.72 


2850 


-47.55 


2066 


-18.95 




3600 


-21.47 


2950 


-47.34 


2070 


-19.51 




3 700 


-18.51 


3075 


-46 . 64 


207 4 


-19.31 




3800 


-16.43 


308 8 


-46 . 81 


2078 


-19.03 




3900 


-20.92 


3188 


-47 . 17 


2082 


-20.04 




3975 


-18.47 


3325 


-47.14 


2086 


-20.06 




4125 


-21.14 


3362 


-47.27 


2090 


-19.76 




4225 


-26.53 


3500 


-47 .24 


2094 


-19.12 




4325 


-24.27 


3663 


-46 . 85 


2098 


- 19.03 




4438 


-25.00 


3763 


-46.63 


2102 


-19.01 


- 


4525 


-21 . 41 


3838 


-46 .00 


2106 


-18.64 




4625 


-20.66 


3938 


-45.49 


2110 


- 19.25 




4713 


-23. 11 


4088 


-45.19 


2114 


- 19.41 




4813 


-21.31 


4150 


-45.55 


2118 


- 19.48 




4900 


-22 . 86 


4188 


-46.10 


2122 


- 19.56 




5025 


-29. 27 


4288 


-46 . 43 


2126 


-18.28 




5162 


-32.61 


4500 


-46.21 


2130 


-17.17 




5262 


-28 . 85 


4588 


-46.76 


2134 


-17.76 




5338 


-28 . 89 


4662 


-47.82 


2138 


-19.66 




5400 


-29.68 


4750 


-48 . 27 


2142 


-21.11 




5500 


-27.05 


4913 


-48 . 34 


2146 


19.75 




0 


-22.19 


5000 


-48.96 


2150 


19.75 




0 


-19.87 


5075 


-48 . 89 


2154 


19.75 




0 


-19.65 


5188 


-48 . 86 


2158 


19.75 




0 


-20.04 


5238 


-49.20 


2162 


19.75 




0 


-23.01 


5363 


-49.70 



55 



NPS 



UPPER HYDROPHONE 



ANALYIS METHOD! NORMAL ANALYSIS 

DATE OP ANALYSIS! 12/1/86 



RECORDED MAGNITUDE LEVELS (dB re IV) 



DELAY( MS ) 


DIRECT MAGN 


DIRECT PRQ( FTZ ) 


REFT* MAGN 


REEL FRQ( HZ ) 


NOISE MAGN 


1900 


19.75 


0 


19.75 


0 


-47.29 


1916 


19.75 


0 


19.75 


0 


-49.07 


1932 


-9.14 


775 


19.75 


0 


-46 . 88 


1948 


-10.98 


638 


19.75 


0 


-47 . 79 


1964 


-13 . 94 


787 


19.75 


0 


-47 . 37 


1980 


-11.66 


800 


19.75 


0 


-46.11 


1996 


-9.63 


963 


19.75 


0 


-46 . 98 


2012 


19.75 


0 


19.75 


0 


-46 . 98 


2028 


19.75 


0 


19.75 


0 


-46 . 98 


2044 


19.75 


0 


19.75 


0 


-50.59 


2060 


-7.81 


1388 


-13.54 


725 


-41.32 


2076 


-10.37 


1738 


-15.72 


1025 


-40. 88 


2092 


19.75 


0 


19.75 


0 


-40.88 


2108 


-10.00 


1738 


-15.14 


1063 


-40.72 


2124 


-12.50 


2050 


-12.39 


1325 


-38.68 


2140 


19.75 


0 


19.75 


0 


-38.68 


2156 


-12.03 


2088 


-9.50 


1375 


-39.19 


2172 


-13.90 


2350 


-13 . 87 


1713 


-41.26 


2188 


19.75 


0 


19.75 


0 


-41.26 


2204 


-13.13 


2325 


-14.15 


1688 


-41.00 


2220 


-16 . 85 


2613 


-13.45 


1913 


-39.81 


2236 


-14.03 


2438 


-14.45 


1775 


-41.69 


2252 


-15 . 76 


2787 


-11.38 


1975 


-39.08 


2268 


-13.21 


2975 


-19.75 


2225 


-42.37 


2284 


19.75 


0 


19.75 


0 


-42.37 


2300 


19.75 


0 


19.75 


0 


-42. 37 


23 16 


19.75 


0 


19.75 


0 


-42.37 


2332 


-15 . 87 


3275 


-22.67 


2588 


-40.30 


2348 


-17.11 


3500 


-20.52 


2825 


-42.85 


2364 


-15.06 


3362 


-20.24 


2662 


-41.32 


2380 


-17.99 


3650 


-17.01 


2988 


-39.05 


2396 


19.75 


0 


19.75 


0 


-39.05 


2412 


-16.95 


3688 


-20.10 


3000 


-41.72 


2428 


-18.05 


3950 


-18.58 


3225 


-38.69 


2444 


19.75 


0 


19.75 


0 


-38.69 


2460 


-17.19 


3975 


-18.88 


3250 


-41.44 


2476 


-18.27 


4263 


-29.36 


3513 


-39. 46 


2492 


19.75 


0 


19.75 


0 


-39.46 


2508 


-17.26 


4288 


-23.68 


3650 


-41.00 


2524 


-25.22 


3925 


-25.22 


3925 


-40.76 


2540 


19.75 


0 


19.75 


0 


-40.76 


2556 


-17.02 


4600 


-21.30 


3900 


-39.97 


2572 


-18.99 


4888 


-19.34 


4188 


-38.58 


2588 


-17.62 


4688 


-21.09 


4012 


-39.92 


2604 


-18.86 


4963 


-20.25 


4238 


-39.87 


2620 


-15.51 


5262 


-22.99 


4538 


-40.13 


2636 


-18 . 87 


5012 


-24.54 


4288 


-40.12 


2652 


-15.35 


5288 


-22.95 


4563 


-40.37 


2668 


-17.57 


5500 


-22.17 


4850 


-41.76 


2684 


19.75 


0 


19.75 


0 


-41.76 


2700 


-16.64 


5475 


-24.00 


4725 


-41.09 


2716 


19.75 


0 


-22.76 


5062 


• -40.32 


2732 


-24.83 


575 


-25 . 79 


5225 


-43 . 41 



56 



NPS 



MIDDLE HYDROPHONE 



ANAL Y I S METHOD: NORMAL ANALYSIS 

DATE OF ANALYSIS: 12/1/86 



RECORDED MAGNITUDE LEVELS ( dB re IV) 



DELAY( MS ) 


DIRECT MAGN 


DIRECT PRQ(HZ) 


REEL MAGN 


REFL FRQ(HZ) 


NOISE MAGN 


1930 


-7.18 


500 


19.75 


0 


-43.05 


1946 


-5.68 


750 


19.75 


0 


-40.55 


1962 


-7.43 


963 


19.75 


0 


-37 . 96 


1978 


-6.07 


787 


19.75 


0 


-42 . 38 


1994 


-6.92 


1100 


19.75 


0 


-37.34 


2010 


- 2.66 


1313 


-13.44 


750 


-34.46 


2026 


19.75 


0 


19.75 


0 


-34.46 


7042 


-4.24 


1388 


-11.47 


825 


-35.67 


2058 


-11.41 


1713 


-7.71 


1275 


-35.05 


2074 


-8.44 


1400 


-12.00 


912 


-35 . 40 


2090 


-9.51 


1775 


-10.34 


1213 


-33 . 98 


2106 


19.75 


0 


19.75 


0 


-33.98 


2122 


- 12.61 


1862 


-10.84 


1313 


-36.46 


2138 


-15 . 45 


2150 


-14.66 


1575 


-39.82 


2154 


-12.28 


1988 


-7.44 


1375 


-35 . 86 


2170 


-16.65 


2238 


-14.10 


1688 


-36.95 


2186 


-12.23 


2138 


-12 . 53 


1475 


-38.29 


2202 


-20.14 


2450 


-12.91 


1775 


-36 . 94 


2218 


- 14.61 


2737 


-23.38 


2150 


-38.87 


2234 


- 16.73 


2550 


-14.09 


1925 


-38. 19 


2250 


- 11.35 


2800 


-19. 81 


2113 


-35.76 


2766 


-13.70 


2963 


-30.08 


2400 


-40.76 


2282 


-11.45 


2787 


-18 . 13 


2088 


-39.75 


2298 


-14.29 


3050 


-30.94 


2463 


-34.28 


2314 


- 14.28 


3325 


-24.50 


2725 


-39.68 


7 330 


-13.08 


3050 


-26.08 


2425 


-40.02 


2346 


-13.65 


3337 


-27.21 


2712 


-37 . 96 


2362 


-15.95 


3563 


-21 . 88 


3013 


-37.62 


2378 


-12.60 


3387 


-22.94 


2800 


-39.30 


2394 


-15.89 


3738 


-25.15 


3125 


-36 . 58 


2410 


-18 . 46 


3925 


-18.28 


3463 


-36.68 


2426 


- 15.35 


3775 


- 21.60 


3113 


-36.10 


2442 


- 15.15 


4037 


- 17 . 45 


3387 


-36.16 


2458 


- 14.81 


3850 


-16.75 


3312 


-38. 12 


2474 


- 14.23 


4150 


-20.62 


3538 


-36.54 


2490 


-13.49 


4425 


- 17.55 


3825 


-36.87 


2506 


- 13.19 


4238 


-18.34 


3650 


-37.63 


2522 


- 13.81 


' 4475 


- 21.49 


3950 


-35.05 


2538 


■ 14.54 


4725 


- 19.42 


4150 


-37.45 


2554 


• 12.97 


4538 


- 17.91 


3925 


-36.71 


2570 


- 14.72 


4813 


- 19.53 


4250 


-34.68 


2586 


-15.36 


5062 


- 20.28 


4413 


-38. 15 


2602 


13.83 


4850 


-17.80 


4300 


-36 . 48 


2618 


- 14.22 


5188 


- 17.22 


4550 


-34.21 


2634 


- 14.23 


5012 


-17.36 


4388 


-36.44 


2650 


-13.65 


5238 


-18. 16 


4650 


-35.15 


2666 


-17.83 


5500 


-22.87 


4975 


-39.41 


2682 


-15.08 


5262 


-18.82 


4700 


-37 . 40 


2698 


-18.71 


5500 


-22.69 


5012 


-36.37 


2714 


-47.10 


3300 


-22.79 


5212 


-40.03 


2730 


- 46 . 87 


4163 


-23.76 


5138 


-40.44 


2746 


-45.35 


3288 


-24.53 


5350 


-37.55 


2762 


19.75 


0 


-19.77 


5188 


-39.64 



57 



NPST 



LOWER HYDROPHONE 



ANALYIS METHOO: NORMAL ANALYSIS 

DATE OP ANALYSIS i 11/30/86 



RECORDED MAGNITUDE LEVELS ( dB re IV) 



DETAY( MS ) 


DTRECT MAGN 


DTRECT FRQ(HZ) 


REEL MAGN 


REPL PRQ(HZ) 


NOTSE MAGN 


1950 


- 14 . 98 




662 


19.75 


0 


-51.50 


1966 


-14.68 




912 


19.75 


0 


-48.51 


1982 


-11.24 




1188 


19.75 


0 


-48.03 


1998 


-13.12 




975 


19.75 


0 


-50.42 


2014 


“9.93 




1263 


-21.03 


575 


-47.79 


2030 


“ 14.56 




1388 


-18.47 


813 


-47.43 


2046 


-9.72 




1-925 


-21.52 


675 


-47.14 


2062 


-14.58 




1563 


-16.53 


875 


-46.68 


207 8 


- 18.08 




1775 


-19.52 


1250 


-44.77 


2094 


“14.41 




1700 


-20.90 


1050 


-46.75 


2110 


- 15.23 




1900 


-15.40 


1313 


-45.23 


2128 


19.75 




0 


19.75 


0 


-45 . 23 


2142 


-16 . 86 




2075 


-13.21 


1388 


-44.71 


2158 


19.75 




0 


19.75 


0 


-44.71 


2174 


-17.68 




2213 


-17.51 


1475 


-46.16 


2190 


-21.85 




2525 


-15.29 


1813 


-46. 43 


2206 


-18.56 




2288 


-22.28 


1700 


-44.55 


2222 


- 19.50 




2625 


-20.51 


1913 


-46.75 


2238 


19.75 




0 


19.75 


0 


-44.94 


2254 


-17.98 




2725 


-20. 10 


1975 


-46. 54 


2270 


-16.99 




2963 


-24.49 


2263 


-47.05 


2286 


-16.51 




2825 


-23 . 31 


2075 


-46. 45 


2302 


-16.05 




3125 


-25.68 


2475 


-46.14 


2318 


16 . 32 




3050 


-27.78 


2400 


-46. 23 


2334 


-17.49 




3350 


-26.05 


2712 


-45 . 83 


2350 


-18 . 22 




3250 


-29.23 


2613 


-45 . 74 


2366 


-23 . 30 




3588 


-23 . 90 


3025 


-46.55 


2382 


-18.01 




3450 


-26.33 


2813 


-45 .79 


2398 


-20.50 




3763 


-18.76 


3075 


-45.33 


2414 


-20.27 




3650 


-22 . 28 


3013 


-45.67 


2430 


-21.14 




3963 


-19.63 


3225 


-44.09 


2446 


-20.07 




3838 


-21.37 


3113 


-45 . 34 


2462 


-20.42 




4188 


-27.08 


3438 


-44.97 


2478 


-21 . 12 




406 3 


-23 . 38 


3350 


-45.16 


2494 


- 21.38 




4388 


-28.26 


3788 


-44.19 


2510 


19.75 




0 


19.75 


0 


-44. 19 


2526 


-20.11 




4488 


-25.96 


3825 


-44.24 


2542 


-io.77 




4775 


-25. 18 


4050 


-42.58 


2558 


- 19.76 




4575 


-25.24 


3888 


-43.69 


2574 


- 19.69 




4900 


-25 . 83 


4200 


-43 . 90 


2590 


“20.20 




4750 


-26.31 


4125 


-♦.22 


2606 


-18.95 




5100 


-27.98 


4388 


— 43 . 


2622 


-20.06 




5000 


-24.44 


4338 


-42 . 88 
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